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We develop here a procedure to obtain regular static configurations as resnlting from dynamical 
gravitational collapse of a massive matter cloud in general relativity. Under certain general physical 
assumptions for the collapsing cloud, we find the class of dynamical models that lead to an equilib¬ 
rium configuration. To illustrate this, we provide a class of perfect fluid collapse models that lead to 
a static constant density object as limit. We suggest that similar models might possibly constitute 
the basis for the description of formation of compact objects in nature. 
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I. INTRODUCTION 

We know today that in nature compact objects such as white dwarfs and neutron stars do form from gravitational 
collapse of massive stars that finally settle to a regular, stationary configuration. Typically, such a compact source 
that results from collapse is supported by nuclear or electromagnetic forces in terms of internal pressures, but we 
do know that gravity plays a crucial role both in the process of collapse itself as well as in the final equilibrium 
configuration. This is true also for larger entities in the universe, such as galaxies and clusters of galaxies that may 
have a very massive compact core at the center. Therefore it is useful to investigate the role played by gravity in the 
formation of compact gravitating objects. The general theory of relativity constitutes the natural framework for such 
a study whenever the gravitational field becomes strong (i.e. when the size of the object approaches its Schwarzschild 
radius). 

In fact, static regular perfect fluid compact sources in general relativity have been used to model very dense objects 
such as neutron stars for a long time now (see for example [l| or and references therein). Furthermore a theory of 
relativistic elasticity provides the first step towards a physically meaningful description of the inner structure of such 
objects 0. Also relativistic models for compact sources of several kinds have been proposed as description of yet to 
be found astrophysical exotic objects. These vary from boson stars (see 0 and references therein), to gravastars (see 
i): preon stars (see for example @) and quark stars (see for example 0)- Even though their existence is not proven 
the search for theoretical models that describe exotic compact objects is very important since it is related to many of 
the unsolved problems of modern astrophysics, such as dark matter, dark energy, the formation of structures in the 
universe and the final fate of gravitational collapse of massive stars. 

While the astrophysical aspects of such general relativistic models are very much under discussion and further 
investigation is needed, we certainly know that the physical process that gives rise to neutron stars in nature is the 
gravitational collapse of a massive star, somewhere in the range of 4-40 solar masses. In such a collapse process 
that takes place at the very end of the life-cycle of a massive star, the outer layers of the star are blasted away in 
a supernova explosion while the inner heavy core implodes to form a neutron star, a very high density object of the 
radius of the order of 10 kilometers. 

As mentioned above, there have been attempts to model such static massive compact objects within a relativistic 
framework, both analytically and via numerical simulations. Analytical models are particularly interestiM as they 
provide insights on the properties that the interior geometry of such objects must have (see for example [8|). Also, 
from the study of such models it is possible to derive important information on the emergence of negative pressures 
and violation of energy conditions that are expected to happen close to the core when the boundary of the object 
approaches or surpasses the event horizon. This provides some insights on the nature of the equation of state for 
extremely dense matter purely from a relativistic perspective (see for example 0). However, due to the intrinsic 
difficulty of Einstein’s field equations, the general relativistic description of dynamical collapse itself that settles to 
such a final static configuration has not been explored so far in detail. Such a process, and especially the late stages 
of the gravitational collapse are the phases where the general relativistic effects should certainly be important, and 
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in fact these could be the dominant ones which really rule the final outcome of collapse. It is therefore very useful to 
investigate, within the framework of the gravitation theory, the scenarios where a dynamical gravitational collapse of 
a massive matter cloud would evolve to result into an eventual regular static configuration. In this direction in M 
it was shown that large equilibrium static configurations can be obtained from gravitational collapse, when certain 
pressures are present. In [ij the observational features of such exotic compact objects were studied. Such studies are 
very important from the point of view of astrophysics as they can tell us how to distinguish an exotic compact object 
from a black hole. For example, in the case of gravastars it has been shown that a realistic model must necessarily 
have anisotropic pressures [l2| and that the spectrum of quasinormal modes of such objects is considerably different 
from that of a black hole of the same mass |1^ . 

In the present work we would like to ask and consider the question if one can obtain general relativistic collapse 
that, neglecting the presence of other forces, still settles to a static regular configuration. We find the answer 
to be in the affirmative, and we present here a class of perfect fluid collapse models tending to an equilibrium 
configuration which is described by a well-known static interior. The perfect fluid collapse model here, although 
lacking a constitutive equation relating pressure and density, can be considered physically viable since it satisfies 
all basic physical requirements (such as energy conditions, regularity, continuity equation for the matter fields) and 
develops from a physically realistic initial state to a well known and physically viable final state. 

In Section 2, we give the general procedure to generate an equilibrium configuration as the final state of the 
gravitational collapse of a perfect fluid. Then Section 3 illustrates this procedure with a specific application showing as 
to how the constant density equilibrium objects result from collapse. The last section summarizes concluding remarks 
which try to bring out the relevance of gravitational collapse in general relativity towards generating compact static 
objects in nature. In the following we make use of geometrical units for which G = c = 1. 


II. STATIC CONFIGURATIONS FROM GRAVITATIONAL COLLAPSE 

The most general spherically symmetric spacetime describing a gravitationally collapsing matter cloud takes the 
form. 


ds^ = —-I- -^dr^ + E^dVi^ , (1) 

Gr 

where </>, R and G are functions of the comoving coordinates t and r. For a matter source made of perfect fluid the 
energy-momentum tensor is given by, Tq = —p, = Tjf = p. The density and pressure are then coupled to 

the metric functions via Einstein equations, which in this case can be written as, 


p = 

F 

R^R ’ 

( 2 ) 

p = 

F' 

R^R' ’ 

( 3 ) 

<^' = 

p' 

P + P ’ 

( 4 ) 

G = 

4 ^- 

( 5 ) 


where (') denotes a derivative with respect to r and (') denotes derivatives with respect to t. 

In the above, the Misner-Sharp mass F{t,r) describes the amount of matter enclosed by the shell labeled by r at 
any given time t, and is given by, 

F = R{l-G + e-‘^‘^R^) . ( 6 ) 

We note that since the system of Einstein equations has six unknowns, namely p, p, (j), G, F and R in five equations, 
there is the freedom to choose one free function. Once such a choice is made, for example by supplying the mass 
profile F for the collapsing cloud, then the full collapse evolution and the final state of collapse is determined fully 
by the Einstein equations above. On the other hand, equivalently, providing an equation of state for the collapsing 
matter that relates the pressure to the energy density also closes the system completely. 

In order for the collapse model to be physically realistic we must require a few physical reasonability conditions. 
These include absence of shell crossing singularities (which implies the condition R' >0), the weak energy condition 
that imposes positivity of the energy density (namely p > 0) and of the sum of density and pressure (p -I- p > 0), 
regularity of initial data and also absence of trapped surfaces at the initial time from which the collapse develops. 
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Finally the collapsing cloud must in general be matched to a generalized Vaidya exterior spacetime at the boundary 
Rb{t) = in the case where there is a mass inflow or outflow from the star during the collapse phase M- In 

the case where the pressure vanishes at the boundary and the Misner-Sharp mass within the boundary is conserved, 
the matching can be done with an exterior Schwarzschild spacetime [l^. 

Here we take without any loss of generality, that the dependence of F in f goes through R in the form F = F{r,R), 
which is always possible whenever R is monotonic in t (as is the case for any gravitational collapse). If we want to 
describe collapse the prescription I? < 0 is enough to ensure that the cloud is not expanding and if i? = 0 is reached 
asymptotically, this ensures the desired monotonic behaviour. In this case we can perform a change of coordinates 
from (r, t) to (r, R), thus considering t = t{r, R) (in the following (, r) will be used to denote a derivative with respect 
to r in the {r,R) coordinates, so that F' = + F,rR')- Using the additional gauge freedom that comes from the 

scaling of the model, we set the initial time so that R{r,ti) = r. 

To solve the system of Einstein equations then we must hrst of all fix the behaviour of the free function that is left 
in the system. In the present case we choose this to be the mass function F(r, R). Supplying the mass function then 
fixes the behaviour and the future time evolution of the collapsing matter cloud. Then equations m and m give 
p{r, R) and p{r, R) as functions of R. Integration of equations dH) and dS]) then gives the metric functions (f) and G as: 


R) 

G{r,R) 


P' 


lo P + P 


df, 




(7) 

( 8 ) 


where &(r), which is usually called the velocity profile for the collapsing shells, is another free function coming from 
the integration and which is related to the kinetic energy of the infalling matter shells. The system is then solved 
once we integrate the Misner-Sharp mass equation m, written in the form 


t.R = - 


G-l 


to obtain the function t{r,R), or equivalently the physical radius of the cloud R{r,t). 
With this, the metric describing the collapsing spacetime is then given as. 


= —e p+p‘^^dF + 


R'- 


b{r)e^ 7 , 


® ^dR 

r R' 


dF FR^dW. 


(9) 


( 10 ) 


It is clear of course that in general it might not be possible to fully integrate the system of Einstein equations, but 
that may not always be needed also. 

Our main aim here is to construct the dynamical collapse of a massive cloud that settles to an equilibrium state in 
which the pressure balances the gravitational attraction. In order to achieve this we must choose the mass function 
and the velocity profile suitably. In fact in general the dynamics as implied by the Einstein equations can lead to 
three possible final outcomes for a continual gravitational collapse: Firstly, an indefinite complete collapse (as in the 
case of a dust cloud, for example), where all matter falls into the central spacetime singularity, which may be covered 
in an event horizon or may possibly be visible to an external faraway observer. Secondly, there may be a bouncing 
behaviour where the infalling matter shells re-expand after reaching a minimum radius. Finally, the collapse may 
settle to obtain a static final object, in which case we must balance the velocities and pressure of the matter content 
of the cloud in order to obtain the limiting behaviour that lies between the two earlier outcomes. 

The equation of motion dSD can be written as an effective potential for any fixed r as. 


Vir, R) = -if 




( 11 ) 


If i? < 0 at all times and F = 0 is not reached even asymptotically, the collapse then does not halt and all the matter 
is crushed into the final central singularity resulting in the formation of a black hole or a naked singularity [l^ . If the 
final outcome is to be a static, regular configuration then some conditions must be imposed on the pressure profile so 
that the collapse will eventually halt, reaching zero velocity and acceleration. The conditions that must be imposed 
so that the metric evolves towards an equilibrium configuration therefore are, 

R = R = 0, (12) 


and they are equivalent to U = U/j = 0. 
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It is clear of course that a static configuration cannot be achieved for the dust i.e. pressureless collapse models, where 
V is negative at all times. Nevertheless from the analysis of the potential when non-zero pressures are introduced, we 
see that 14, as a function of R for any fixed shell r, can have in general two zeros. Then the static configurations can 
be obtained from those potentials for which both the zeroes coincide. This implies that V has a local maximum for 
which V — 0. By linearizing the potential close to the equilibrium point we can see that the static configurations can 
be achieved in the present comoving coordinate system asymptotically as t grows to infinity. 

If the solution of the equation of motion ([9]), given by i?(r, t), tends asymptotically to an equilibrium solution Re{r) 
such that the conditions given by equations m are satisfied, then the collapse evolves towards a static equilibrium 
configuration. 

In order to obtain the equilibrium scenario we must then choose the free function F (r, R) during collapse in such 
a way that the quantities F, (p, and G tend to their equilibrium limits, namely F{r,R) —?> Feir) = F{r,Re(r)) and 
similarly for p and G, where (peif) is given by the Einstein equations and Ge{r) is determined by the imposition of 
the static conditions m- 


Ge{r) = 1 - ^ . 

(G.fl)e = GAr,Re{r)) = ^-^ 


(13) 

(14) 


Note that the velocity profile b{r) appearing in equation Q has been absorbed into Ge{r). The above equations 
are direct consequence of the equilibrium conditions (1121) . The key point here is the following: Given any static 
equilibrium configuration, which is characterized by specifying the mass function Fe(r) of the static system, the free 
function F{t,r) during the collapse is to be chosen in such a way that F{t,r) Fe{r) in the limit. The class of all 
such dynamical collapses as specified by this condition will then settle to the equilibrium limit. 

At the equilibrium we can define a new radial coordinate from Re{r) = z and rewrite the metric functions as. 


Fe{r) = F{z), pe(r) = p{z), Ge{r) 


GC) = l-f 


Then from equations m and we recover the Einstein equations for a static interior, given by. 


(15) 


pi.z) = ^, p,z = -{p + p)(t>,z, (16) 

where the second equation is the well-known Tolman-Oppenheimer-Volkoff equation. Then the third static Einstein 
equation, namely 


P = 




(17) 


can be obtained from equation ([2]) once we impose the equilibrium conditions and we make use of equation (O at 
equilibrium. 

The spacetime geometry given by the metric (1101) during collapse, once the equilibrium is reached, reduces to the 
familiar static spherically symmetric geometry given by 


Ct 




(18) 


which can be matched to a Schwarzschild vacuum exterior at the boundary Zb = Re{rb)- We note that the equilibrium 
conhguration again is fully determined by the choice of one free function, which we take to be the mass function F{z). 
An important point to notice here is that the above metric in principle need not be necessarily regular at the center 
since any singularity that might eventually form is obtained as the result of continual collapse from a regular initial 
data. 

Generally, when studying static perfect fluid sources of the Schwarzschild spacetime, one requires a set of physically 
viable conditions to be satisfied [l7|. These are the usual energy conditions, the matching conditions with exterior 
Schwarzschild geometry given by the vanishing of the pressure at the boundary, a monotonically decreasing behaviour 
for the energy density and pressure, and regularity at the center. As we have seen, this last requirement of regularity 
at center in the static final state could in principle be omitted. Finally one wishes to impose also that the sound speed 
in the cloud be smaller than the speed of light, thus requiring p/p < 1. 

We note that in general if energy conditions are satisfied during collapse they will be satished by the equilibrium 
configuration as well, namely the positivity of the energy density and sum of density and pressure at the origin follows 
from the same condition during collapse. Further we note that requiring only the weak energy condition allows for 
the possibility of some kind of negative pressures. 
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III. CONSTANT DENSITY INTERIORS AS COLLAPSE LIMIT 


As an example of the procedure described above to generate regular static equilibrium configurations resulting as 
final state of a dynamical collapse, we now consider a well-known static interior for the Schwarzschild metric which is 
given by a constant density distribution. This is the first interior solution obtained by Schwarzschild himself and it 
is simple enough to illustrate how a static source can be achieved from a collapsing scenario as a limit. Other static 
interiors, such as those obtained by Tolman [l^ . can be investigated in the similar manner but we will not go into a 
detailed analysis for other spacetimes here. 

The main question we address here is: Can we obtain the above constant density regular Schwarzschild interior 
equilibrium configuration as limit of a regular and physically viable dynamically collapsing matter cloud, within the 
framework of general relativity? 

In order to have constant density for the interior metric we take 

F{z) = , (19) 

such that 

p{z) = pq . (20) 


The system of static Einstein equations can easily be integrated and the pressure is then given by 


P(z) = Po 


y/T^- yl-C 


i-c^-Sv/r^ 


( 21 ) 


where C = ^ is given by the boundary value for F{zb) = 2M such that it matches a Schwarzschild manifold with 
mass parameter M at p{zb) = 0. 

From the above condition we see that p(0) = and we must impose 3^1 — C < 1 (corresponding to the 

upper mass bound known as Buchdahl limit M < 0 ) in order for the central pressure to be finite and positive 

[^ . Furthermore if we require that the sound speed within the cloud, defined by = p/p, does not exceed the speed 
of light we must require a further constraint on C. For example, evaluating the sound speed at the center we get the 
condition C < 3/2 (corresponding to a further upper mass bound M < jZb). 

The static metric m then becomes the well-known constant density Schwarzschild interior given by, 

-3\/r^| dt^+ ^ Fz^dn'^. (22) 

1-C4 

Now we look for the classes of dynamical collapse models that lead to the above static metric as the limit of collapse. 

All those dynamical collapses with F{r,R) going to an equilibrium configuration limit Fe{r) = ^Re{r)^ and 
satisfying the equilibrium conditions m and (imi will asymptotically tend to the above static spacetime. This is in 
general a wide class of spacetimes, and not all the possible choices of F will give a physically valid dynamical collapse 
evolution. Of course, in general it will not be possible to fully integrate Einstein equations during collapse to give all 
the global dynamical solutions tending to the static limit. 

Therefore, of the class of all possible F for collapse we now choose one specific class of models where p does not 
vary in time, thus choosing p(r, R) = p{r) = Pe{r). The reason to do so resides in the fact that with this choice we are 
able to integrate equation (I^, thus obtaining an exact analytical solution for collapse almost fully, while in general 
exact integration of Einstein field equations for collapse proves to be unattainable. From equation (ED then we must 
have. 



F{r, R) = y{r)R^ + w{r), 


(23) 


where the freedom to specify F is reflected in the freedom to choose y and w, and the pressure is then given by 

p(r) = -3y{r). (24) 

In such a case, as the cloud shrinks from the initial radius R{r,ti) to its final configuration Re{r), the radial profile 
for pressure as a function of the comoving radius r remains unchanged. From equation m we therefore must take 


y{r) 


./TTc- 

o^_ _ y _ H 

3 - 3\/r^ 


(25) 
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and by requiring that the energy density goes to a constant value po in the limit of equilibrium we must choose, 


w{r) 


2po 




l-C 


Rejr)^ 


-3^/T^ 


:Re{rf. 


(26) 


This finally fixes the choice of the free function during collapse. The boundary of the cloud (chosen such that 
Rei^b) = Zb where the pressure vanishes) corresponds to the physical radius i?(rb, t) that shrinks from the initial value 
Ri^bRi) = Tb to the equilibrium value Zb- Note that since p{rb) = 0 during collapse, we must have F{rb, Rb{t)) = 
w(rb) = 2M which means that the cloud must match to a Schwarzschild exterior. The energy density during collapse 
is then given by 


P = 


rR^ 


■ w, 


R^R' 


-P- 


(27) 


The above choice of the pressure profile during collapse was made so that equation ([5]) can be integrated. Performing 
the integration we obtain the metric function G as. 


G{r,R)=b{r) 


f y,rr^ + W,r 


2 


(28) 


where b{r) is a free function coming from the integration and we have considered here the initial scaling condition to 
be R{r,ti) = r. Therefore by imposing the equilibrium conditions we must take the free function 6(r) to be 

, . , _ Rejr) - yRejrY - W / y^rRejrY + 

Re{r) V y,rr^+W^r 

This fixes the dynamical evolution completely. It is easy to see that all quantities obtained depend upon the unknown 
function R{r,t), which still needs to be determined. To obtain the full solution for the dynamical collapsing cloud 
one should in principle integrate the two remaining equations, namely ([1]) and ([S]). Equation can be written as. 





3y,rR^R' 

y^rR?' + W,r ’ 


(30) 


and its integration gives the other metric function 4>{r,t) also in terms of R{r,t). The only equation that finally 
remains to be solved then is the equation of motion (jS]) which can be written as. 


R = -e^ 


+ G-1. 


R 


(31) 


We can treat this equation as an ODE for the function R of t for each fixed radial shell r. We then note that this 
is a first order ODE in t that, given the continuity of the functions involved, always admits a solution. Therefore 
integration of this equation gives the metric function i?(r, t) and thus solves the system of Einstein equations entirely. 
The dynamical spacetime m tending to the static limit given by the metric (1^^ is then hnally written as. 


ds^ = —e 


2 / 0 ’' 


3y,: 


■ r^r' 


rR'^ + v. 


dr 


dv 


R'‘^dR 


b(r) ( 


y ,rR'^,1 


+ R^dn^ 


(32) 


Equation (l27ll can be used to study the behaviour of the central density during collapse. If we suppose that all the 
functions involved are finite and continuous in the radial direction, we can then expand the area coordinate R{r, t) 
near the center as 


R{r,t) = ai{t)r + a 2 {t)r'^ + ■■■ (33) 

with the initial conditions ai{ti) = 1 and aj{ti) = 0 for j > 1. Then as collapse evolves R{r,t) will tend to the 
limit R{r,t) — > Reij") with Re{r) = aieV + a 2 ef‘^ + ••• from which we can see that the central pressure will have the 
following behaviour 


= « + (34, 

From the above equation for the central density we obtain a series of constraints that can be imposed on the model 
in order for it to be physically valid: 
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• Positivity of initial density: 


2\/r^ 


(1-aL) 


> - 1 . 


• Initial density smaller than final density: 




(l-aL)<0. 


• Initial sound speed smaller than I: 


I - - 2al^VT^ 


< 1 . 


Note that the last constraint above implies that imposing the speed of sound within the cloud to be smaller than 1 
at the initial time is sufficient for the same to be smaller than I during the whole evolution. 

Further to this, one would like to require that there are no shell-crossings during the collapse evolution. This means 
that the shell labeled by ri and r 2 do not overlap at any point during the evolution. The singularities that originate 
due to shell crossings (namely when R' — 0) signal a breakdown of the coordinates at that point and therefore imply 
the impossibility to describe the future evolution of the cloud using that coordinate system. These are generally 
believed to be weak singularities that have no physical relevance being due only to the coordinate choice |2l| . The 
requirement that there are no shell-crossing singularities is given by i?' > 0 which, by making use of equations ® 
and equation can be written as. 


R' = 


{p + p)R'^ 


> 0 . 


(35) 


It follows then that provided that the weak energy condition is satisfied throughout collapse, the condition for 
avoidance of shell-crossing is just given by > 0. Then it is easy to check for the above toy model that if 
the condition is satisfied at equilibrium it will be satisfied throughout collapse. This may be considered as a very 
reasonable physical condition in that it is basically the requirement that the mass should be increasing with the 
increasing radial coordinate r. 

We note that imposing during collapse physically realistic conditions such as regularity of initial data, avoidance of 
shell crossing, sound speed smaller than the speed of light, allows us to restrict the array of possible scenarios that can 
evolve from collapse. This eventually allows us to narrow down the large number of existing static interior solutions 
on the basis of those that can be achieved via realistic collapse. 


IV. CONCLUDING REMARKS 

It is clear for many years now, that the physics of gravity and the general theory of relativity must play a crucial 
role when it comes to understanding very compact objects in nature. This has been the key motivation for the many 
attempts over past years for a fully general relativistic modeling and description for such objects. These models also 
greatly help in understanding the accretion disk properties and related observational details for such objects and the 
associated high energy phenomena, with observations now offering an increasing degree of preciseness and detail. 

Since the original interior solution with constant density was provided by Schwarzschild in 1915, there have been 
a lot of studies on matter sources in general relativity. Static interiors sustained by perfect fluid pressures were also 
investigated by Tolman [Tsil who found a number of solutions describing compact objects. Interiors sustained by 
anisotropic pressures were investigated in and [2^, while interiors sustained only by tangential pressures were 
firstly studied in [^ . Since then the zoo of static interiors to the Schwarzschild vacuum spacetime has become very 
crowded (see for an overview of interior solutions with perfect fluids). In principle, given a generating function 
in the form of (j)^z, it is possible to construct any interior matching smoothly to the vacuum Schwarzschild metric 
[ 2 ^. Nevertheless, as pointed out by Lake and Delgaty [ 2 ^, most known solutions fail to fulfill some of the physical 
requirements listed above. Furthermore, such solutions need not arise naturally from a physically realistic collapsing 
conhguration and their connection with collapsing models has not been thoroughly investigated. From the point of 
view of astrophysics the critical densities that lead to the formation of different kinds of compact objects were first 
studied in and have been thoroughly investigated ever since. Nevertheless our knowledge of the states of matter 













that are possible in extreme conditions is still very limited and it is possible that, together with the neutron star 
matter equation of state, other kinds of ‘exotic’ fluids might exist in nature, possibly for more compact and dense 
objects. For this reason many theoretical models for ‘exotic’ compact stars have been proposed over the years. What 
we really need to decide is which ones will be the most relevant models to choose from, from a physical perspective, 
out of such a wide collection and variety offered by general relativity. 

Under the situation, while a general relativistic description of such static compact objects where strong gravity 
effects are highly important is not only essential but is also inevitable, what is really necessary is to isolate the models 
which would be physically most relevant and significant and the equations of general relativity must necessarily be 
the starting point in this search. It is in this spirit that we argued here that the static configurations that arise 
from a general relativistic gravitational collapse, which we already know to be a very relevant physical process in 
the universe, could provide us with a physically preferred set of models, as opposed to the plethora of very many 
other static solutions to Einstein equations as mentioned above. For example, an important question that needs to be 
addressed is whether the gravastar models described in [l^ can be obtained via a collapse mechanism such as the one 
described here. In such a case negative pressures and non-ideal fluid models are likely to play a significant role as well. 
Also modifications to general relativity in the strong field may be required to balance the gravitational attraction. A 
simple and generic toy model for relativistic collapse leading to the formation of such exotic compact objects is still 
missing. Such a model would play a role similar to that of Oppenheimer-Snyder collapse for black hole physics [^ . 

In other words, we propose here that the compact static configurations as arising from gravitational collapse, as 
per the procedure outlined here, may play a physically more significant and essential role, and that these should be 
investigated in further detail. Further work in that direction in under way. We showed here how it is possible to 
obtain perfect fluid static interiors from dynamical collapse with regular initial data. This shows that analytic models 
that have been considered to describe the equation of state for neutron stars, such as the Tolman interiors (see for 
example Hi), can be obtained quite naturally from gravitational collapse. Further physical considerations are needed 
if one wishes to apply similar models in realistic astrophysical scenarios, and the study of different profiles for densities 
and pressures will also be important in this connection. 
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